etf->Draw("pow(86.329723/(tof1),0.8)*sqrt(dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)-2.8*(86.329723/(tof1)))/((86.329723/(tof1))*(86.329723/(tof1))-log((86.329723/(tof1))/(sqrt(1-(86.329723/(tof1))*(86.329723/(tof1)))))))*69.438048-1.2552140:1.606*sqrt(1.-pow(86.329723/(tof1),2))/(86.329723/(tof1))*(1.+(pdcTaUWire[1][0]-12.84*x_f1)/8500.)*pow(20/z0*1.606*sqrt(1.-pow(86.329723/(tof1),2))/(86.329723/(tof1))*(1.+(pdcTaUWire[1][0]-12.84*x_f1)/8500.)/2., 1/60.)-0.044>>(1000, 1.95, 2.3, 1000, 14, 39)", "pdcTaUWire[1][1]!=-999.&&pdcTaUWire[1][0]!=-999.&&x_f1!=-999.", "col")



et->Draw("33.3786*(sqrt(dE[1]*exp(0.00035*(pdcTaDWire[1][1]+40.))))-18.2236:beta>>(500, 0.45, 0.62, 500, 14, 39)", "pdcTaDWire[1][1]!=-999.&&beta>0.", "col")


etf->Draw("beta0*sqrt(dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.))/(beta0*beta0-log(bg0)))*exp(-1.*beta0):aoz0>>(1000, 1.95, 2.25, 1000, 0.25, 0.73)", "pdcTaUWire[1][1]!=-999.&&bg0>0.", "col")




etf->Draw("sqrt(dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.))):tof1>>h1(500, 135., 165., 1000, 0.4, 1.4)", "pdcTaUWire[1][1]!=-999.&&pdcTaUWire[1][0]!=-999.&&x_f1!=-999.", "col")



86.329723/(tof1)*
sqrt(
  dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)) /
    ((86.329723/(tof1))*(86.329723/(tof1))
      - log((86.329723/(tof1))/sqrt(1-(86.329723/(tof1))*(86.329723/(tof1)))))
)

// z v.s. tof1
etf->Draw("sqrt((dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)))/(0.00228*pow(tof1-121.43353,2)+6.973)):tof1>>(1000, 135, 165, 1000, 0.23, 0.63)", "pdcTaUWire[1][1]!=-999.", "col");
// z v.s. aoz
etf->Draw("sqrt((dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)))/(0.00228*pow(tof1-121.43353,2)+6.973)):aoz0>>(1000, 1.95, 2.25, 1000, 0.23, 0.63)", "pdcTaUWire[1][1]!=-999.", "col")
// z v.s. tof1 -> local
etf->Draw("sqrt((dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)))/(0.00228*pow(tof1-121.43353,2)+6.973)):tof1>>(300, 139, 147, 300, 0.5, 0.62)", "pdcTaUWire[1][1]!=-999.", "col")
// tof corrected for Kr
etf->Draw("sqrt((dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)))/(0.00228*pow(tof1-121.43353,2)+6.973))-(tof1*tof1*0.000830073-0.238031*tof1+17.06):tof1>>(300, 139, 147, 300, 0.5, 0.62)", "pdcTaUWire[1][1]!=-999.", "col")

etf->Draw("sqrt((dE[0]*exp(0.0015*(pdcTaUWire[1][1]+30.)))/(0.00228*pow(tof1-121.43353,2)+6.973))-(tof1*tof1*0.000830073-0.238031*tof1+17.06):aoz0>>(1000, 1.95, 2.25, 1000, 0.23, 0.63)", "pdcTaUWire[1][1]!=-999.", "col")